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ABSTRACT 


This study is concerned with the problem of 
constructing an optimal military deployment plan for 
sealift assets during a period of conflict. The deployment 
problem is formulated as a set-partitioning optimization 
problem with a minimax objective. An algorithm for solving 
this problem is presented and it is based on solving a 
sequence of related, but simpler, linear programming 
problems by the column generation technique. The results 
of the model are ship schedules to meet the cargo 
requirements of the deployment plan in aminimum amount of 
time. Various implementation strategies are discussed as 
well as the occurrence of integer solutions. In addition, 
computational experiments for several small to medium size 


examples are presented. 
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I. INTRODUCTION 


A. PROBLEM STATEMENT 

During a military deployment, troops, equipments and 
Supplies must be transported from ports of embarkation to 
ports of disembarkation. Generally, the standard modes of 
transportation used in this operation are trucks, trains, 
airplanes and_= ships. Because of the limited amount of 
available resources and transport assets, planning becomes 
essential for a successful deployment. During peacetime, 
cargo transportation can be routinely scheduled and the 
normal criterion for a deployment plan is its cost (or 
operating expense). However, during a period of conflict 
(or crisis), expenses become secondary and it 1S more 
important to transport the troops and cargoes to their 
destinations as fast as possible. 

This study restricts itself to the problem of 
constructing an optimal deployment plan which employs only 
sealift assets. Many of today's naval deployment plans are 
constructed manually and in an ad hoc manner. This process 
1s quite time-consuming and does not guarantee to produce 


even a near optimal plan. 


B. BACKGROUND 

Research in deployment planning for both industrial and 
military applications has been concentrated on constructing 
deployment plans which minimize operating costs. In a 
survey article, Ronen [Ref. l:pp. 119-126] describes the 
various modes of operation for cargo ships and provides a 
classification scheme for ship routing and scheduling 
models. In amore recent article, Brown et al. [Ref. 2:pp. 
335-346] present and solve the crude oil tanker scheduling 
problem formulated as an elastic set-partitioning model. 

On the military side, Goodman [Ref. 3] formulates the 
problem of scheduling the naval surface combatants of the 
Atlantic Fleet as a generalized set-partitioning problem. 
The resulting constraint matrices in both formulations of 
Brown et al. [{Ref. 2] and Goodman [Ref. 3] have a large 
number of columns which must be generated beforehand and 
correspond to all feasible ship schedules. In a Naval 
Postgraduate School master's thesis, Collier [Ref. 4] 
formulates the deployment planning problem employing four 
standard modes of transportation as a linear programming 
problem, and solves it by the MPS III Mathematical 
Programming System developed by Ketron Management Science, 
Inc. [{Ref. 5]. Related to Collier's study, Lally [Ref. 6] 
uses the General Algebraic Modelling System, GAMS [Ref. 7], 
to solve the problem of minimizing the number of sealift 


assets required to carry out a given deployment plan. 


C. OBJECTIVE 

In previous formulations of the deployment or ship 
scheduling problem, the primary objective is to minimize 
cost which is the most appropriate for peace-time military 
operations and for industry. This thesis addresses the 
same problem, but with a different objective: to minimize 
the duration of the deployment. Ines Darticular,  1t 
considers the construction of schedules for sealift assets 
to transport cargoes from their ports of origin to their 


ports of destination in the minimum length of time. 


II. PROBLEM FORMULATION 


To formulate the crisis deployment problem, the 
following data are assumed to be given: 


1. The ports of embarkation and disembarkation for each 
cargo 


2. The distances between ports 
3. The number of ships with their speed 


4. The compatibility between each ship and each cargo 


When a ship 1S compatible with a cargo, we mean that 
the ship is compatible with both the cargo and its ports of 
embarkation and disembarkation. Therefore, in constructing 
the compatibility information one has to consider, for 
example, the ship and cargo type as well as the ship draft 
and the channel depth of both ports. 

It 1s assumed that all cargoes are configured into full 
shiploads. This implies that when a ship picks up a given 
cargo, it must deliver it before any other cargo can be 
picked up. Therefore, the ship must travel to the port of 


disembarkation directly from the port of embarkation. 


Ae. MATHEMATICAL FORMULATION 
The problem of scheduling sealift assets ina crisis 
Situation can be formulated as a variation of the standard 


set-partitioning model as follows. 


immarces: 


= indexes shiploads of cargoes, where 
i=1,2,...,m and m is the number of 
shiploads; 

J oc indexes ships, where j=1,2,...,n and n 


is the number of ships available; 


k - indexes a ship schedule. 
Data: 
Sjk —- a binary vector representing the kth feasible 
schedule for ship j. The it® component of Sjx, 


denoted by Sijk, 18S: 


1, if the kt® feasible schedule 


for ship j includes shipload i; 
Sijk = < 


0, otherwise; 


tyk - the completion time of schedule Sjx. 


(The calculation of tjyx is described in a 


numerical example below.) 


Decision variables: 


1, aif the kt® feasible schedule 
for ship j 1s selected for the 


Xjk = < deployment; 


O, otherwise. 


Problem P11: 


min [{ max ( & tik *Kik 7s. 0) fp euike eee ) j 
k k 


subject to: 


n 
> = Siyk *® Xjk a | for i=]. om (15 
j=1 -k 
x Xk $s 1 for j=l 4a=.,n (2) 
k 
Xjxn € fO7 2)e 
The term cK tyk ®*Xjk in the objective function 


represents the time for ship j3 to complete its assigned 
schedule. Therefore, max [{ fx tik*Xik,...,2k tux *Xox | 
represents the completion time of the longest schedule in 
the deployment plan. Since the deployment 1s considered 
completed only when all cargoes are delivered to their 
destinations, the completion time of the longest schedule 
becomes the length of the deployment plan, which 1s to be 
Minimized. The first set of constraints in Problem Pl 
ensures that all cargoes are picked up by at least one ship 
and the second guarantees that at most one schedule is 
assigned to each ship. 

In addition, the objective function of Problem Pl is a 


nonlinear convex function since it is a point-wise maximum 


ot a set of linear functions. However, Problem Pl can be 


transformed into a linear problem as follows. 


Problem P2: 


Minimize 2Z 


subject to: 


nN 
>» = Sijk * Xju ae b Hor 1 = ei: 
j=1 &k 
x Xk <1 for j=l son: 
k 
Cas yo Gk 2 0 ie) ce 3 | eee lS 
k 


Siiicee (Oe lse 


Note that the last set of constraints defines the objective 


mumection of Problem Pl. 


B. AN EXAMPLE 

To illustrate the above integer programming 
formulation, consider the deployment problem depicted in 
Bugure 2.1. There are two ships, Ship 1 and Ship 2, 
available for the deployment. Initially, Ship 1 is docked 
at Port Pland Ship 2 at Port P2. There are 3 shiploads 
whose ports of embarkation and disembarkation are given in 


Table 2.1. Assume that both ships are compatible with all 


ports and cargoes. The lines connecting various ports in 
Figure 2.1 represent possible movements between ports. The 
TABLE 2.1 


A LIST OF POE AND POD FOR THE SHIPLOADS 


SHIPLOAD 





3 
ike = P] 
OT KS ee 
Je 


PII 
2 may CC 


—- TRAVEL TIME 


(In days) 


Figure 2.1 Data for the Deployment Problem 


numbers adjacent to the lines represent the travel times 
for both ships, l1.e., they have the same speed (assumed 
constant regardless of cargo loading). 

Consider a schedule for Ship 1 which includes picking 
up cargoes 1 and 3. The binary vector representing this 
schedule has 3 components (Since there are 3 shiploads) 
with the first and the third components having the value 
one and the second component the value zero. To carry out 
this schedule, Ship 1 can use one of the two possible 
routes: one which picks up cargo 1 first and then cargo 3 
and the other which picks up the cargoes in the reverse 
order. Using the time given in Figure 2.1, the first route 
requires 8 days to complete and the second requires only 7 
days. Since the objective 1S to minimize the completion 
time of the longest schedule, the completion time of 7 days 
1s assigned to this schedule, i.e., tia = 7. In general, 
the completion time tjx is the time for ship j to carry out 
schedule k using the shortest route. Tables 2.2 and 2.3 
display all possible schedules along with their completion 
times for Ships 1 and 2, respectively. Note that the 
schedule discussed above is the schedule Si4 in Table 2.2. 

The optimal deployment plan for this example consists 
of two schedules: Si2 for Ship 1 and S24 for Ship 2, and 
requires 7 days to complete. In terms of decision 
variables, xi2 and x24 equal one and all other variables 


equal zero. 


The explicit formulation of the above example is given 


in Figquremz:. 2. 


TABLE 2.2 


POSSIBLE SCHEDULES FOR SHIP 1 











SCHEDULES 


poe | se | oe 





COMPLETION 
TIME 
(IN DAYS) 


suman 7 | ef 
mmowa [epee pepape pa 
wmoo> fe fet=[=fe*{=]s 







TABLE 2.3 


POSSIBLE SCHEDULES FOR SHIP 2 












SCHEDULES 


Ses | Ser | See | Sm | 
TIME eZ 


2 f) 
(IN DAYS) 


moose fefefs[sfels 
mes ste fs fefet= ts {a 
moo ste fefst=tet= |. 






COMPLETION 







10 


minimize 2 

subject to: 

Xii t+ Xi4+tXist Xi7 + X21+ X24+X25+ x27 
X12+ XistXi6etXi17 + X22+ X25 +X26+X27 


X13 +Xi4+ Xie+Xi7 + X23 +X24+ X26 +X27 


X11 4X1 2 +X1 34X14 4+X154+X16+X17 


X21 X22 X23 ¢X24 4X25 tX26 +X27 


z- ( tir*Xritti2*xXi2ttis *Xigttia*Xiattis*Xisttis *Xi6etti7*xi7 | 


z- { tei *xX2itte2*xXe2ttes *Xeattea *X2attes *X25tte6*x26tt27*xX27 ] 


X11, X12, X13, X14, X18, X16, Xi7 € [0,1] 


A21, X22, X23, X24, X28, X26, X27 € {0,1] 





Figure 2.2 A Formulation Of The Example Problem 


al 


III. A SEQUENTIAL SOLUTION TECHNIQUE 


The solution procedure presented below addresses 
Problem P1 (and P2) indirectly. This procedure takes 
advantage of the fact that there exists a simpler problem 
which is related to Problem Pil (and P2). By sequentially 
solving a number of these simpler problems, one can arrive 


at a solution to Problem Pi (and P2). 


A. A RELATED PROBLEM 

In certain situations, it is not so critical that the 
planner obtains a deployment with the minimum duration. 
Instead, the duration of the deployment, say t days, has 
been set by the top command and the planner only has to 
find a feasible plan which can be completed within this 


given length of time. To formulate this problem, define: 


K; (t) = { k : Sjx is a feasible schedule for ship j 


and tae << =a 


That is, Ky (t) is the set of schedules for ship j which 
can be completed within tt days. Then, we have the 
following problem which we refer to as the feasibility- 


seeking problem. 


WZ 


Problem P3(t): 


min Wi 


uns 


af 


subject to: 


n 
wi + X ( x Siu * Sabee) 2 1 POreal fe 2—i, o..«,m? 
j=1 keKy (tT) 
>» Xj k = ip forma. j=l,<6.5n;: 
keKy (tT) 
Wier, hime 10,1). 


where wi 1S an auxiliary variable to indicate whether or 
not shipload i will be left undelivered by the deployment 
plan. If the optimal solution to Problem P3(t) is greater 
than zero, i1t means that t is infeasible. In this case, 
one or more. shiploads must be left undelivered or 
additional assets are required to obtain a plan which can 
be completed int days or less. Thus, if tis a feasible 
duration, Problem P3(t) will produce a feasible plan. 

Note that Problem P3(t) is parameterized by t. If the 
Minimum duration for a deployment plan, t*, is known, then 
the solution to Problem P3(t*) solves Problem Pl (and P2) 
as well. Otherwise, by varying t and resolving Problem 
P3(t) ln a systematic manner, one can obtain a solution to 
Problem Pl (and P2). A strategy for searching for the 
Minimum feasible duration t* is discussed in Chapter IV. 


aN 


To illustrate the feasibility-seeking problem, consider 
the deployment problem presented in Chapter II. Assume 
that the planner is told to construct a plan witha 


duration of 8 days. 


Then, 
Ki(8) = ({ 13025053 ceo and 
Kz (8) = { 2 2% 3G, 4, om 6 ie 


that is, the eligible schedules for this plan with a 


completion time of 8 days or less are those listed Tables 


3.1 and 3.2. In this case, the optimal objective function 
value for Problem P3(8) is zero, because 8 days is a 
feasible duration. Each of the following pairs of 
schedules’ for Ships 1 and 2: (Si1,S26), (Si12;)Sz240m 
(Sis ,S25 7, (Sia ce and (Sie ,S21), constitutes a 


deployment plan that can be completed within 8 days. 
Similarly, if one solves Problem P3(t) with t equal to 

7 days (the optimal duration), the optimal objective 

function value is still zero, and the pairs (Si2,S24) and 


(Sia,S22) are the only feasible deployment plans. 


B. A COLUMN GENERATION APPROACH TO THE FEASIBILITY-SEEKING 
PROBLEM 

Since the feasibility-seeking problem searches for a 

feasible deployment plan and does not have a real objective 

function to optimize, one expects that the relaxation of 

the integrality restriction would not produce fractional 
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sSomutions too often. This observation is corroborated by 
the computational result presented in Chapter IV in which 
integer solutions are obtained for over 90 per cent of the 
problems. Henceforth, we treat Problem P3(t) as a linear 
programming problem. 

As a linear program, Problem P3(t) has many columns. 
To avoid generating these columns a priori, we apply the 
column generation technique, 1.e., the Dantzig-Wolfe 
decomposition method [see, e.g., Ref. 8], to Problem P3(t), 


and the following decomposed system is obtained. 


TABLE 3.1 


ELIGIBLE SCHEDULES FOR SHIP 1 WHEN t=8 





SCHEDULES 






COMPLETION 
eM Gs 
(IN DAYS) 


a 
moos [ets fete fs 
moos fete fs[=[s 






1S 


1 Pes er 


ELIGIBLE SCHEDULES FOR SHIP 2 WHEN t=8 









SCHEDULES 


COMPLETION 
hie 
(IN DAYS) 


moos [= fefef= |e {e 
mmoos{e f= fete|sis 
moose fe[=fste[s 







Master Problem (MP): 


min 


Wi 


wm 3 


a 


SUD JeECE sto: 


n 
oy 3 Sigk * Xjyk + Wi 2k for i=l... . meee 
j=1 kebyj (tk) 
D2 Xj k <2 for j=1)...,neae 
keLy (tT) 
0 Ss xi Ss for all a7] 
Subpreoblem foreship. 3 (2SP1 (a) aoe 
m 
kK =Aerarg min { vi + 3% Si gieee* Soe 
keKy (tT) i1=1 


where ui is the dual variable corresponding to constraint 
set (1), i.e., the cargo (shipload) constraints, and v; is 
the dual variable corresponding to the constraint set (2), 
1.e., the ship constraints. We refer to ui: as the it 
cargo dual and vj as the jt® ship dual. 

The column generation technique starts with an initial 
set of feasible schedules, Ljy(t), for each ship j. This 
initial set Ly (t) may be an empty set. The master problem 
is solved andthe dual variables ui and vj are obtained. 
From this set of cargo and ship duals, one or more 
subproblems are solved thereby generating additional 
schedules (columns), Sjx’, which are subsequently added to 
the set Lj (tT). The master problem is then resolved with 
the additional schedules (columns) and the new cargo and 
ship duals’ ars ained. The cycle then continues until 
ene objective t£unction value of Problem SP(3) 1s 
nonnegative for all j, i.e., all schedules have nonnegative 
reduced cost. This signifies that optimality is achieved. 
Figure 3.1 illustrates the cycling between the master and 
subproblem in the column generation technique. 

As stated above, the subproblem is unnecessarily hard. 
In theory, it is not necessary to add schedules (columns) 
with the most negative reduced cost to the master problem. 
Any schedules (columns) with negative reduced cost would 
suffice. The following subproblem produces negative 


reduced cost schedules for the master problem. 
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Master Problem 


m 
minimize = OW 
Steal 


subject to: 


L Sigk * Xjk + Wi 


jJ=1 = keLy (tT) 


L Xjk 
keLy (T) 


O S$ Xjx 


The subproblem 
produces a new 
Columns Sika, 
for the master 
PEODleIN> a..6:.5 
Sjk’ 1s added Subproblem 
Oey Capa ees 


min 


keKkjetw 


for 15) 22.7 





The master 
problem produces 
cargo and ship 
duals for the 


subproblem. 





Figure 3.1 - The Column Generation Technique. 


1s 


Subproblem SP2(j): 


For ship j find an index k' such that k'‘ ¢e Kj (t) and 


s 
+ 
kms 


Sigk *® Ui < 6. 


i=1 


If k' solves Problem SP2(j), k*' is an acceptable 
schedule. For details concerning the generation of 
acceptable schedules, the reader is referred to a related 


Master's Thesis by LCdr Svein Buvik [Ref. 9]. 
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IV. IMPLEMENTATION AND COMPUTATIONAL RESULTS 


To implement the column generation procedure we 
modified the revised simplex code described in Ref. 10, to 
solve the master problem. In this modification, we allow 


the algorithm to restart from the last optimal solution 
after one or more new schedules (columns) have been added 
to the master problem. Since the set partitioning problem 
is usually degenerate, we also reinvert the basis at every 
ten iterations. As for the subproblem, we employ the 
algorithm developed by Buvik [Ref. 9]. Both the master and 
subproblem algorithms are written in FORTRAN 77 and 
compiled by the IBM VS FORTRAN compiler. All runs were 
performed on an IBM 3033 AP computer at the W.R. Church 


Computer Center of the Naval Postgraduate School. 


A. PROBLEM DATA 

For our experimentation below, we consider the 
deployment scenario in which cargoes must be moved from the 
ports along the east coast of United States to ports in 
Europe. Table 4.1 lists approximate distances between 
various ports. The number of shiploads for our deployment 
problems are varied between 5 and 50 and the list of all 50 
shiploads along with their POE's and POD's are given in 
Table 4.2. The number of ships assigned to the deployment 
are assumed to be between 2 and 30 ships and the initial 
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location of all 30 ships are given in Table 4.3. The speed 
of all 30 ships are between 18 and 25 knots, and on the 


average a ship is compatible with 75 % of the shiploads. 


TABLE 4.1 


DISTANCES BETWEEN PORTS 
(IN NAUTICAL MILES) 


rs ee ee 
EL oe 











where 

Ny. = New York 
Norf. = Norfolk 
Charl. = Charleston 
Jax. = Jacksonville 
Ham. = Hamburg 
Wilh. = Wilhemshaven 
Rot. = Rotterdam 
Antw. = Antwerpen 
Chbe = Cherbourg 
Pens. = Pensacola 
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TABLE 4.2 
LIST OF SHIPLOADS 


SHIPLOAD POE POE SHIPLOAD POE POD 
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. EE 
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~ 
oO ~~] uw Wi DO] 
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ay 


a 
ml Ww 
CEEEE 


Od 
O| wow] of a] om 
ele CEREEE +f efofeled = 


to 
ra 
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ROT MT M7] bd 
Ol] i Wy wo 
. 


POE POD 
1. New York 1. Hamburg 
2. Norfolk 2. Wilhelmshaven 
3. Charleston 3. Rotterdam 
4. Jacksonville 4. Antwerpen 
5. Pensacola 5. Cherbourg 


ez 


TABLE 4.3 
INITIAL DISTANCES BETWEEN SHIPS AND PORTS 
(IN NAUTICAL MILES) 


ER 
ase [eas [580 | 720 | tis 
[400 [355 [450 | 620 | 250° 
3280 [945 [650 | ta30 | 290 
SE 
/_350_ [320 | 990 | 3900 | 2400” 
[sf 00 [380 [650730 | 9280 
SC 
[=| 200 [00 | 400 | s00 | a0” 
[#0 | | s00 | 700 | 900 | 2400 
[ax [aso [00 [00 [875 [1078 
[aa [380 [98 [200 [375 | 78 
[a3 [80 | aso [a0 [ies [700 
[as [#80 [00 [00 [265 | 900 
SS 
[30 [720 [avs [aes | ae0 | soo 
[32 an07 | a500_ | s00[ 305 [3300 
[22 [ 920 [ars [365 [200 | 400 
[33 | #50 [350 [300 [400 | 00 
[34 [“4350_|a2s0_[ 200 [3300 [i600 
[as [aas0[ 975 [700 [600 | 2330_ 
[3s | aso [275 [700 | s00 [3100 
(a7 #80 [ 675 [00 [300 | 300 
[—35_| 900 [70 | s00 | 400 | S00 
[30 | a0 | a75 [soo [00 | 100” 
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B. STRATEGIES FOR GENERATING SCHEDULES 

As described in Chapter III, the decomposition process 
iteratively solves the master and subproblem in sequence. 
After having just solved the master problem, the subproblem 
obtains the cargo and ship duals from which it generates 
one or more negative reduced cost columns. At this point, 
there are several possibilities regarding the ship(s) for 
which the subproblem should generate schedules (or 
columns). The first obvious strategy is to generate a 
schedule for Ship 1 in the first iteration, a schedule of 
Ship 2 in the second iteration and so on until a schedule 
for each ship has been generated. At which point, the 
cycle of generating schedules (columns) begins again with 
Ship 1. The second strategy 1s to generate schedules for 
ship in the descending order of the ship duals, and the 
third strategy is just the reverse of the second strategy, 
i.e., generates schedules in the ascending order of the 
ship duals. The other strategy, which has been considered 
and soon after discarded, generates schedules for all ships 
during each iteration. This strategy tends to generate the 
same schedule for all ships, which seems redundant since no 
two ships can have the same schedule at optimality. In 
fact, there must exists a solution in which no two ships 
are assigned to pick up the same shipload. Based on this 
observation and preliminary experiments, the last strategy 


is discarded. 
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To compare the three strategies discussed above, we 
solved three feasibility-seeking problems at various 
lengths of deployment, t. The first problem has 30 
shiploads and 20 ships, the second problem has 35 shiploads 
and 20 ships, and the last has 40 shiploads and 20 ships. 
In Figure 4.1, we plotted the average cpu times on these 
three problems against the length of deployment. The first 


and third strategies clearly dominate the second. 


C. SOLVING THE MINIMAX PROBLEM 

As mentioned in Chapter III, one can solve the minimax 
problem by sequentially solving the feasibility-seeking 
problem in the following manner. First, pick an initial 
value for +t and then solve the feasibility-seeking problem 
at this value t. If the optimal objective function 1s 
zero, then the value of t is decreased and the feasibility- 
seeking problem is resolved at the new value. Otherwise, 
the optimal objective function is positive, the value for T 
is increased and the feasibility-seeking problem is 
resolved. 

The efficiency of the above algorithm is clearly a 
function of the initial value for +t. If the initial value 
for t 1s close to the optimal, the feasibility-seeking 
problem has to be solved less often. Thus, it 1S important 
that a good initial value for +t is used to start the 


process of increasing and decreasing the value tl. 
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Figure 4.1 A Comparison of Three Strategies for 


Generating Schedules 
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One lower bound estimate is given by the following 


equation: 
tT. = integer part of [(2*ntr - 1) * trmin + itmin] /spmax] 
where 
ntr = average number of cargoes (shiploads) per ship, 
es. 
number of cargoes 
0 ot i — — ti la le ' 


number of ships 


trmin = the minimum travel distance between POD's, 


itmin = the minimum distance between ships' initial 
positions and POE's, and 


spmax = maximum speed among all ships. 


To understand this bound, note that for each shipload 
assigned to aeship, the ship has to first deliver the 
shipload to its destination and return to pick up the next 
shipload on its schedule. Therefore, the ship has to make 
two trips (or ocean crossings) back and forth between POE's 
and POD's for each shipload, except for the last one for 
which the ship only has to make one trip from a POE toa 
POD. Thus, if ntr shiploads are assigned to one ship, it 
has to make (2*ntr - 1) trips. Since the minimum distance 
between a POE and a POD is trmin, the minimum total 
distance traveled by each ship is (2*ntr - 1) * trmin + 


itmin. The first term represents the distance for trips 
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between POE's and POD's and the second term represents the 
distance from ship's initial position to the first POE. 
Then, dividing the total by the maximum speed among the 
ships gives a lower bound for the optimal t. Table 4.4 
displays the value of the lower bound estimate and the 
correspond values of 1t*, the optimal duration, for 35 
problems. On the average, tT. underestimates t* by 40 &. 

If historical data, e.g., data from previous deployment 
exercises, are available, the lower bound estimate tT can 
be improved by using linear regression. For example, using 
the data from Table 4.4, we obtain the following linear 


equation 


Test = 15.57 + 0.8 * wt 


where Test represent the linear estimate of t* based on T1.. 
Figure 4.2 displays the linear estimate of t* graphically. 
Since linear regression minimizes the squared error, some 
Test naturally overestimates t*. Based on T. and Test, we 
implemented the following search algorithm for tT. 

In the algorithm below, the initial estimate, 1, of 
the optimal duration, t™*, is obtained by taking a convex 
combination of the lower bound and the linear regression 
estimates. It is assumed that the convex weight a, is 
chosen so that 1%: underestimates t*. (Note that this is 
always possible by letting a equals one.) The parameter 6 


equals one time unit which is one day in all our examples. 
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TABLE 4.4 


LIST THE LOWER BOUND ESTIMATES AND ACTUAL VALUES 


ee) 
dled. 33 


OF OPTIMAL T 


PROBLEM NUMBER OF NUMBER 
NUMBER SHIPLOADS OF SHIPS 
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Figure 4.2 





TEST = 15.57 + lore 
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ct 
Lower Bound Estimate versus Actual Values of 


Optimal t. 
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Algorithm 


Step 0: Set ae meee tee Ce Tent and 
set k=l. 

Step 1: Solve the feasibility-seeking problem, Problem 
P3(tk), by the column generation technique. 

Step 2: If the optimal objective function value equals 
O, stop; tr 1S optimal. Otherwise, set 


T+1 = Tr + 6 and k = kK +1. +Go to step 1. 


In Step 2, the current estimate, tT, of the optimal 
duration 1S increased by amount 6. In this manner, the 
current estimate tx approaches the optimal duration t* from 
below and all of the previously generated schedules remain 
feasible to the feasibility-seeking problems in the 
succeeding iterations. One topic for future research is to 
relax the assumption that 171 must underestimate t* and 
allow Tk to be adjusted in either upward or downward 
Girection in Step 2. Table 4.5 summarizes the 
computational results for the above algorithms. [In all 


cases, the value for ais 0.7. 
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TABLE 4.5 


COMPUTATIONAL RESULTS FOR THE COLUMN GENERATION TECHNIQUE 
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D. PERCENTAGE OF INTEGER SOLUTIONS 

Since the approach taken in solving the crisis 
deployment problem is the linear relaxation of the minimax 
set-partitioning problem, it is of interest to investigate 
the question concerning the integrality of the obtained 
results. In theory, the linear relaxation of the problem 
does not always produce an integer solution, in which case 
an integer programming algorithm such as the branch and 
bound method must be employed. However, Table 4.6 
demonstrates that an integer programming algorithm is 


necessary for less than seven percent of the problems. 


TABLE 4.6 


NUMBER OF PROBLEMS WHICH TERMINATE WITH OPTIMAL 
INTEGER SOLUTIONS 


# OF PROBLEMS| # OF PROBLEMS WITH 
15 SOLVED INTEGER SOLUTIONS 


ico e” = 1.29% t* 


Menor 8-8). 1o* Tt 
morc  — 1.09*t* 
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V. CONCLUSIONS AND FUTURE STUDIES 


This study formulates a crisis deployment problem as a 
set-partitioning problem with a minimax objective. An 
algorithm is developed for solving this problem. The idea 
underlying this algorithm is to solve the minimax set- 
partitioning by solving a sequence of simpler, but related, 
feasibility-seeking problems. Each time the feasibility- 
seeking problem produces a better solution for the minimax 
problem. The feasibility-seeking problem is similar in 


form to the minimax problem and both have a large number of 


columns. So to solve the feasibility-seeking problem, the 
column generation technique (as in the Dantzig-Wolfe 
decomposition method) is employed. The computational 


results in Chapter IV verify that this method is effective. 

An important by-product of the above development is 
that the feasibility-seeking problem can also answer the 
question: Can all cargoes be deployed to their final 
destinations in t days? A negative answer to this question 
leads to two natural follow-up questions which provide 
interesting areas for future studies: 


1. How many additional ships are required to deploy all 
cargoes in t days? 


2. If no additional ship is available, which cargoes 
must be left behind? 
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Besides the above areas and the one mentioned in 


Chapter IV, the following areas are also worth studying. 


i 


The scenario considered in this study assume that the 
deployment is completed in one phase. In an extended 
period of conflicts, one may want deployment plans in 
several phases (waves). 


Several embellishments to the current model are also 
possible. 


a. Allow the cargoes to arrive at the ports within 
time windows. The current model assumes that all 
cargoes are alwaysS available for transport. 


b. Allow cargoes in partial shiploads and in 
compatibility among cargoes, e.g., ammunition 
should not be loaded on same ship with fuel. 


es Allow Eos nondeterministic delays ng) the 


completion times. These delays are due to 
unfavorable weather and/or enemy blockade. 
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APPENDIX 


FORTRAN PROGRAM 


KKK KKKEKKEKEKKKEKKKEKEEKKEKEKKEKEKKEKEKEKKEKKKKKKKKEKKKKKKKKKKKKKKKKKKKE 


x x 
* * 
x oh  — x 
= PROGRAM DEPLOY = * 
k eeu Se he = eee x 
« : 
* * 
* * 
a Dakice = Zoey Oe aa oboe * 
" x 
* * 
* « 
- Key Variables: % 
* * 
* M - number of constraints; ‘ 
* N - number of variables; = 
* A - real matrix of dimension M by N containing * 
* the coefficients of the M constraints; * 
* B - real vector of length M containing the right * 
* hand sides of the constraints; * 
* C - real vector of length N containing the * 
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coefficients of the objective function; 
XB - basic variables; 
BINV - matrix of dimension M by M corresponding to 


the inverse of the basic matrix; 


IB - set of indices corresponding to the basic 
variables; 

U - dual variables; 

XTIME - duration of the schedule; 

SEQ - sequence of cargoes to pick up; 

ELL - index of the variable leaving the basis; 

K - index of the variable entering the basis; 

SB - search direction; 

SIGB - maximum feasible step size; 

SHIP - ship number; 

OBJ - objective function value; 

MR - movement requirements; 

eT - travel distances from current ship ports to 


POES (ports of embarkation) ; 
TR - travel distances between POEs and PODs 


(ports of disembarkation) ; 


SPD - ship speed; 
TAU - number of days to complete the deployment; 
COMPAT - matrix of dimensions M by M that contains in- 


formation about the compatibility ship-cargo. 


Oy. 


Subroutines 


The subroutines and their objectives are: 


a oLME I 


a eon ee 


eS Ve 


ees Pee. 


S- Salt eR 


ai RINVR SE 


=(RDAY=s 


= oUB ER 


eek. LM 


= RESULT 


- solves 


- step l 


- step 


- step 


2 


5 


the revised simplex method; 
of the revised simplex method; 
of the revised simplex method; 


of the revised simplex method; 


- updates the "B" inverse matrix; 


- inverts the B matrix; 


estimates an initial value for the number 


of days to complete the deployment; 


Key parameters 


generates feasible schedules; 
computes travel times. 


writes the output. 


NLOA - number of full shiploads ctecarcaqces. 
NSH = number of Ships: 
NPOE - number of ports of embarkation; 
NPOD - number of ports of disembarkation; 
Output : The output provides the following information: 


- objective function value, 


- number of simplex iterations, 


- optimal 


(minimum) number of days to complete the 


deployment, 


- optimal primal solution, 
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- optimal dual solution, 
- ships' schedules, 


- sequence of cargos to pick up per ship, and 


- schedules‘ durations. 
Input / Output devices 

Disk (MOVREQ DATA) input device 07 
Disk (TRAVEL DATA) input device 08 
Disk (FSTDST DATA) input device 09 
Disk (COMPAT DATA) input device 11 
Disk (SPD DATA) input device 12 
Disk (DEPLOUT DATA) output - device 10 


KRKKKEKKKKEKKKEKKKEKEKEEKEKEKEKKEKEKEKKEKEKEKKEKEKKKKKEKEKKKKEKKKKKEKKEKKKEKEK 


MAIN PROGRAM 


IMPLICIT DOUBLE PRECISION (A-H,O-Z), INTEGER (I-N) 


PARAMETER( MM = 100, NN = 2000, KK = 2 ) 
DIMENSION B(MM),C(NN) ,SB(MM) ,U(MM) ,BINV(MM,MM) ,IB(MM), 
&WORK (MM) ,XB(MM) , XCOL (MM) ,MR(100,KK) , SPD (MM) 


REA OeeuMhM NN) eR (15,15), 107T(30,15) 
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INTEGER MR,SHIP,XTIME(NN) ,SEQ(NN,MM) , TAU, TAUL,TAUEST, 
&LSEQ (NN) ,NLOA,NSH,NPOE,NPOD,CTSHIP 

LOGICAL COMPAT (MM, MM) 

CHARACTER*13 MOVREQ, TRAVEL, FSTTIM 


COMMON /UNITS/ NIN, NOUT 


& 


C Initialize variables. 


S 
DATA NLOA,NSH,NPOE,NPOD / 3, 2, 5, 5 / 
DATA A,B,C,XB,BINV / 200000*0. ,100*1.D0,2000-G 
«“100*120D0, 10000*0. ODO W 
DATA SB,WORK,U,IB /100*0.D0,100*0.D0,100*0.D0, 1007GF 
DATA SEQ, LSEQ ,SPD / 200000*0,2000*0,100*0.0D0 / 
DATA MR,TR,IT,XTIME /200*0,225%*0.,450*0.,2000*0 / 
NIN = 2 
NOUT = 6 
JOUT = 0 
TAU = 1 
TAUL = 1 
TAUEST = 1 
c 


C Read the data from the data files. 


C 
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Peep Oy ge MR tt, J) 7 0=1,2),i=1,NLOA) 

Pee Oca ein Ld) ut, NOD) ,Il=1 ,NPOE} 
READ (Om (rt (ig) J=), NPOE) ,I=1,NSH) 
READ (11,*) ((COMPAT(I,J) ,J=1,NSH) , I=1,NLOA) 


READ (12,*) (SPD(I) , I=1,NSH) 


C 


C Estimate an initial value for TAU. 


C 
CALL RDAYS (NLOA,NPOE,NPOD,NSH,TAU,TAUL, TAUEST, SPD, 
SR, LT} 
MD = TAU 
WRITE (NOUT, 8000) NLOA,NSH,NPOE,NPOD, TAU 
WRITE (NOUT, 8030) TAUL, TAUEST 
e 


C Convert input to number of columns (M) and number of 


C rows (N) in the "A" matrix. 


M = NLOA + NSH 


ITER = 0 


41 


5000>-. DO 10 1 — ain 
Bit) Seo De 
SB(I) = 0.0D0 
Ut Tja= 0. ODO 
IB(I) = 0 
WORK(I) = 0.0DO 
XB(1) = 1. 0D0 
XCOL(I) = 0.0D0 
DO 20 J = 1,NN 
C(J) = 0.0D0 
XTIME (3) = 0 
A(I,J) = QO. 
SEQ(J,I) = 0 
SEO (do). 0 
20 CONTINUE 
DO 30 K=1,MM 
BINV (27.5) = ©. DO 
30 CONTINUE 
bh8. CONTINUE 
CTSHIE ="0 
SHIP = 0 
K = 0 
N = 2*NLOA + NSH 


@ 


C Generate input for the revised simplex method. 
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Dome 1 = 1, 
Ieee) s Se 8 


DOestae = 1,M 


hese. 0) DINV( tJ) = 1.D0 
a0 CONTINUE 
40 CONTINUE 


DO 60 I = 1,M-NSH 
Ci) =e. OD0 


60 CONTINUE 


C 


C Generate artificial variables. 


2 
DO, 7 Ogs= Lam 
DOTSO K = 1M 
Phd eibe@ee K) A(J,K) = 1. 
80 CONTINUE 
70 CONTINUE 
c 


C Generate surplus variables. 
C 
DO 90 J = 1,M 


DO 100 K = Mt+1,2*M-NSH 
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100 


90 


110 


NORE, 


IF(( J .EQO. (K-M) ) . AND. (J = bE eins 
& A(J,K)= -1. 
CONTINUE 


CONTINUE 


SUM = 0.DO 
DO 110 I=1,M 
U(l) =) =e4n} 
SUM-=°SUM +7C€(IB (1) ) *xei st) 
CONTINUE 
OBJ = SUM 


O CONTINUE 


Strategy to choose for which ship the next schedule 


will be generated 


SHES gon UP 


Lb (SHIP... EOo NSH +e -) VSHLE sa 


Generate columns as needed by the master problem. 


yo 


CALL SUBPR(U,XCOL,TAU,M,N,NLOA,NPOE,NPOD,NSH,MR,TR, 


&IT,A,COMPAT,IB,XB,SHIP,XTIME,K,SPD, SEG) fsck omemontees 


DO 120) = 
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730 


il 26) 


2 


@: 


C 


140 


SUM = 0.0DO 


WO 1305) = 1, M 


SUM = SUM + BINV(I,J) *A(J,K) 


CONTINUE 


SB(I) = SUM 


CONTINUE 


Perform the revised simplex method. 


Cone ote LAA, B,C,48B, BINV,SB,U;WORK,IB,OBJ;N,M,JOUT, 


&K, ITER) 
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(OBJ .LT. 10.0D-4) THEN 
IF(TAU .EQ. MD) THEN 


MeS=9MD + 10 


TAU = TAUL 
GO TO 5000 
END eE 
NT = 1 


DO 140 I = 1,M 
RERCn bie) GT. 2 OD=3 
CONTINUE 
WRITE(NOUT,8010) NT 
IF(NT .EQ. 1) THEN 
COFTre, 1100 


ELSE 
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-AND. XB(T) 


ray Wid bo 


-9 ) NT=0 


IF(ITER .CT. 4000)> Sie: 


GO 703000 
END IF 
END IF 
GO TO 1000 


C 
C Write the results. 
€ 
1100 WRITE(NOUT,8020) TAU 
CALL RESULT(JOUT,XB,U,C,A,1IB,M,N,OBd ,1TBR,scEOn cee 


&XTIME,NSH,NLOA) 


8000 FORMAT(20X,'PROGRAM OUTPUT’ ,/,20X, '==============', 
&//,6X%,12,1X, 'SHIPLOADS',3X,12,1%,  SHiP See eee 


&1X, 'POES',3X,12,1X, *PODS' ,/,6X, ‘INGRAM ee ti eee 


eS VALUE = “PAU I= 9522/7) 
8010 FORMAT(6X,'NT =',12) 
8020 FORMAT(6X,'*** FINAL (OPTIMAL) TAUS—=— 2 
8030 FORMAT(6X,'** TAUL = ‘, 22K,’ , (AURST= ee, 22) 
SLOP 
END 


SUBROUTINE SIMPLX (A,B,C,XB;,BINV;, SB, UpWORK ee, Cou 


&N,M,JOUT,K,ITER) 
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* This subroutine performs the revised simplex method. * 


IMPLICIT DOUBLE PRECISION (A-H,O-Z), INTEGER(I-N) 
PARAMETER( MM = 100, NN = 2000 ) 

DIMENSION XB(MM),B(MM) ,C(NN) ,BINV(MM,MM) ,SB(MM), 
&U (MM) , WORK (MM) , IB (MM) 

REAL A(MM,NN) 

INTEGER ELL,XTIME (NN) 


COMMON /UNITS/ NIN, NOUT 


JOUT = 0 


200 CONTINUE 


PER = TER + 1 


Peesour 8. 1 ) RETURN 


Sri ROTEPS(XB,SB,SIGB,ELL,M,JOUT) 


Peas OUr BO. 2 ) RETURN 


CALL RSTEP3(XB,C,B,BINV,A,WORK,OBJ,IB,ELL,K,N,M, ITER) 
Ore OBIS. LT. L0.0D—4) THEN 

NT = 1 

DO 10 I = 1,M 


PEeb AT) GT. 1. 0D-3 -AND. XB(1I) .LT. .90 ) NT=0 
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10 CONTINUE 
PP (NT EO] ds LHEN 
CTE Rei leRe tas 
RETURN 
END IF 


END IF 


IF( MOD(ITER,10) .EO. 0 ) CALL RINVRT(BINV Bae 


&WORK,M,N) 


CALL RSTEPI(A,C,SB,U,BINV IBV NM KR coum 


GO TO 200 


END 


SUBROUTINE RSTEP1 (A,C,SB,U,BINV)iB Noe er 


* This subroutine performs the step one of the revised 


* simplex method 


IMPLICIT DOUBLE PRECISION (A-H,0O-Z), INTEGER(I-N) 
PARAMETER ( MM = 100, NN = 2000 ) 
DIMENSION C(NN),SB(MM),U(MM) ,BINV(MM,MM) , IB (MM) 


REAL A(MM,NN) 
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20 


108, 


30 


INTEGER XTIME (NN) 


COMMON /UNITS/ NIN, NOUT 


“6 


il 
he 


TOLCON 


JOUT = 0 


Compute the duals. 


DOwsO J=1,M 
SUM = 0.DO0 
IDLO) 2G) Gan 


SUM = SUM + BINV(I,J)*C(IB(TI) ) 


CONTINUE 
eye = — SUM 
K = 0 
VKMIN = 1.D30 


DO 50 I=1,N 


Check if I is in IB. 


DO 30 J=1,M 
trie EO. IB(d) )» GO TO 50 
CONTINUE 


SUM = C(T) 
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DO 40 J=1,M 
SUM = SUM + A(J,I)*U(J) 
40 CONTINUE 
TF (SUM .<GE. VKMIN) GO-TO 750 
VKMIN = SUM 
K = 1 


50 CONTINUE 


TE(VKMIN LE. =-TOLCON) Ge Wemee 
JOUT = i 
RETURN 


60 CONTINUE 


C 


C “BGEm ss oc 


S 
DO 80 I=1,M 
SUM = 0.D0 
DO 70 J=1,M 
SUM = SUM + BINV(I,J)*A(J,K) 
70 CONTINUE 
80 SB(I) = SUM 
RETURN 
END 
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SUBROUTINE RSTEP2 (XB,SB,SIGB,ELL,M,JOUT) 


* This subroutine performs the step two of the revised 


* simplex method 


IMPLICIT DOUBLE PRECISION (A-H,O-Z), INTEGER (I-N) 
PARAMETER( MM = 100, NN = 2000 ) 

DIMENSION XB(MM) ,SB (MM) 

INTEGER ELL 


COMMON /UNITS/ NIN, NOUT 


BPS ee —6 


it 
oO 


ELL 

SIGB = 1.D30 

Be. LOOrl=1 .M 
Poole irene Ss) GO TO 100 
RATIO = XB(I)/SB(I) 
IF(RATIO .GE. SIGB) GO TO 100 
SIGB = RATIO 
ELL = I 

100 CONTINUE 
IF(ELL .EQ. 0) JOUT = 2 
RETURN 


END 


Syl 


SUBROUTINE RSTEP3 (XB,C,B,BINV,A, WORK, Cbd eee 


GN TM, ETRE) 


* This subroutine performs the step three of the revised 


* simplex method 


IMPLICIT DOUBLE PRECISION (A-H,0O-Z), INTEGER(I-N) 
PARAMETER( MM = 100, NN = 2000) 

DIMENSION C(NN) ,XB(MM),B(MM) ,BINV(MM,MM) , WORK (MM) 
REAL A(MM,NN) 

INTEGER ELL,IB (MM) 


COMMON /UNITS/ NIN, NOUT 


DO 10 -1=2,,/M 
10 WORK(I) = A(I,K) 
CALL PHIPRM(BINV, WORK, ELL,M) 
BO 30 1-1 ™ 
SUM = 0.DO 
DOeZO n= 21M 


SUM = SUM + BINV(l oer ty) 


20 CONTINUE 
XB(I) = SUM 
30 CONTINUE 
IB(ELL) = K 
SUM = 0.DO 


Se 


40 


x 


10 


DOn8 0 — 1 


ay 
SUMe= SUMMe C(IB{I))*XB(T) 
CONTINUE 
OBJ = SUM 
RETURN 


END 


SUBROUTINE PHIPRM(BINV,WORK, ELL, M) 


—_e ae eee eee eee eee eee eee eee eee ee ee ee eee eee es eee eee eee ee ee es eee eee eee eee ee ee ee eee eee 


Paice.  OOUDLE PRECISION (A-H,O-Z), INTEGER (I-N) 
PARAMETER( MM =100, NN = 2000 ) 

DIMENSION BINV(MM,MM) , WORK (MM) 

INTEGER ELL 


COMMON /UNITS/ NIN, NOUT 


TOL) (=a ao 
SUM = 0.DO0 
DO 10 I = 1,M 


SUM = SUM + BINV(ELL,I) *WORK(T) 
CONTINUE 
YSUM = DABS (SUM) 


Pec Mernck. TOL) ) GO TO 20 


23 


20 


30 


40 


50 


60 


WRITE (NOUT,8000) SUM 
S LOR 
CONTINUE 
SUM = 1.D0/SUM 
DO-=20 f= 4. 
BINV(ELL,I) = SUM*BINV(ELL,1I) 
IF( (BINV(ELL,I) .LT. TOL) .AND. (Baily Vee 
.GT. -TOL) ) BINV(ELL, 2) )— 70 70he 
CONTINUE 
DO 6040 7=s1.M 
iF (d=. bo. bil) CO. Loree 
TEMES 70 D0 
DOT 40° f= mM 
TEMP = TEMP + BINV(J,1) *WORK(T) 
CONTINUE 
IF( (TEMP .LT. TOL) .AND.. (TEMP <CG2 tele 
TEMP = 0.0D0 
BO 50-15-71 7M 
BINV(J,I) = BINV(J,1) = TEMP*BINV {ERE 
IF( (BINV(J,I) .LT. TOL) .AND. (BINV (aaa 
-GT. -TOL) ) BINV(J,1) = G2Gbe 
CONTINUE 
CONTINUE 


RETURN 
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SoU 0 ehORMAT (6X), 0**** ERROR *** NEW MATRIX WOULD BE’, 
&' SINGULAR, INNER PRODUCT =',G15.6) 


END 


SUBROUTINE RINVRT (BINV,A,IB,WORK,M,N) 


* This subroutine reinverts the basis. 


iMEnle LEmeeOCUBLE PRECISION (A-Hy7O-Z), INTEGER (I-N) 
PARAMETER ( ZERO = 0.0D0, ONE = 1.0D0 ) 

PARAMETER ( MM =100, NN = 2000 ) 

DIMENSION BINV(MM,MM), IB(MM), OMAT(MM,MM), WORK(MM) 
REAL A(MM,NN) 


COMMON /UNITS/ NIN, NOUT 


POL = 1. D-6 
BO-10 If = 1, M 


DOn20 J = 1, M 


BINV(I,J) = ZERO 
OMAT(I,J) = A(I,IB(J)) 
20 CONTINUE 
BINV(I,I) = ONE 
10 CONTINUE 


8, 


C Locate maximum magnitude element on or below the main 
C diagonal. 
‘@ 

DO] 30K = ia 


IF ( K }Oaee eee 


IMAX K 


AMAX 


DABS (OMAT(K,K) ) 
Kel S=sne 
DO 40 I = KP1, M 
IF ( AMAX .LT. DABS (OMAT(1,K))) THEN 
IMAX = I 
AMAX = DABS (OMAT(I,K)) 
ENDIF 
40 CONTINUE 
C 
C Interchange rows IMAX and K if IMAX is not equal to K 
& 
DP CiMAX ONE. K) THEN 
DOl SO e) -= ar 
ATMP = OMAT(IMAX, J) 
OMAT(IMAX, J) = OMAT(K,J) 
OMAT(K,J) = ATMP 
BTMP = BINV(IMAX,J) 
BINV(IMAX,J) = BINV(K,Jd) 


BINV (K,.0 )> =) Bae 
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50 CONTINUE 
ENDIF 


ENDIF 


C 


me Test for singular matrix. 


C 
TEeeCDABS (OMAT(K,K)) sbT. 1.0D-6) THEN 
WRITE (NOUT,8000) K, OMAT(K,R) 
ELSE 
DIV = OMAT(K,RK) 
Be oO J = i, iM 
OMAT(K,J) = OMAT(K,J)/DIV 
ie mCMnAtukee yessbr. TOL) .AND. (OMAT(K,J) 
& 7G le— hOuy ) OMAT(K,3) = 0.0DO 
BINV(K,J) = BINV(K,J)/DIV 
re eeerrny (Kp uy. Gls) TOL) SANDS (BINV(K, J) 
& Pree LOL) BINViK; J) = 0.0D0 
60 CONTINUE 


Dee O le-— 1) M 
AMULT = OMAT(I,K) 
IF( (AMULT .LT. TOL) .AND. (AMULT 
& -GT. -TOL) ) AMULT = 0.0DO 
LES St eeNee Kk) 6oTHEN 
DO 80 J = 1, M 


OMAT(I,J) = OMAT(I,J) - AMULT 


oie, 


& = OMAT (hee 
BINV(I,J) = BINV(I,J) = AMULT 
& x BONY (iqon, 


IF( (BINV(I,d) .LT. TOL) eenire 


& (BINV(I,J3) .GT. -=TOime 
& BINV (1, J) 3 ="O20D0 
80 CONTINUE 
ENDIF 
70 CONTINUE 
ENDIF 
30 CONTINUE 


8000 FORMAT({(* * ERROR: BASIS IS SINGUIAR Pei loec 
RETURN 


END 


SUBROUTINE SUBPR (U,XCOL,TAU,M,N,NLOA,NPOE, NPOD, None 
&MR,TR,IT,A,COMPAT,1IB,XB, SHIP, XTIME, Korey, SEG Gee 


&CTSHIP) 
* This subroutine generates feasible (acceptable) columns * 


IMPLICIT DOUBLE PRECISION (A-H,O-Z), INTEGER {(I-N) 
PARAMETER ( MM = 100, NN = 2000, KK = 2, JJ = 2000 ) 
DIMENSION XCOL(MM) ,U(MM) ,UU(MM) ,V(MM) ,XB(MM) ,SA(MM), 


&IBWMM) , SPD( MM) 


> 


REAL A(MM,NN) ,TR(15,15) ,IT(30,15) 

INTEGER VIND(MM) ,PRED(0:JJ) ,LOAD(O:JJ) ,TIME(0:JJ), 
& MR(100,KK) ,FROLD, TOLD, PATH(0:MM) ,CURLD,COUNT, 
& LENGTH, LASTND, SHIP,TT,MLNGTH,XTIME(NN) ,CTSHIP, 

& CTSON(0:MM) ,LSEQ(NN) ,TAU,STACK(0:JJ) , TOP, 

& SEQ(NN,MM) , CTBACK(0:MM) 

LOGICAL COMPAT (MM, MM) 

DOUBLE PRECISION MIN,MINRC 


COMMON /UNITS/ NIN,NOUT 


c 
C Initialize the variables. 
C 
LIMIT = N 
moOO0 CTSHIP = 0 


2000 DO 10I1I= 1,MM 


SA(I) = 0.0D0 

DUG) == 0. D0 

V(I) = 0.DO0 

VIND(I) = 0 
a0 CONTINUE 


NNEG = NLOA 


MINRC = 0.0D0 


e 


@eeroort the duals. 


3), 


DO 20 I = 1,NLOA 
LSS ie) = Sor) 
20 CONTINUE 
TF( N .GE. (2*NLOA + NSH + 1) eerPee 
DO 30°24 = 714 


DO 40 J = 1,M 


IF( (IB(J) .GT. (2*NLOA + °NSHis eatiie 
& (A((NLOA + SHIP),IB(J)) .NEs92. DO}M—eaAND. 
& (XB(J) 4~GT. 0,5D0) jo TREN 
SA(I) = SA(I) + A(T, IB 
END IF 
40 CONTINUE 
IF( SA(I) .LT. 1.0D0 ) UU(T)> = =2eens 
310 CONTINUE 


END) tere 
C 


C Check ship-cargo and ship-port compatrbita es 


c 
DO 50 I = 1,NLOA 
IF( .NOT. COMPAT(I,SHIP) ) UU(IT) = 99.0D0 
50 CONTINUE 
DO 60 I = 1,NLOA 
MIN = 0.1D-6 
COUNT = 0 
IND = 0 


60 


70 


60 


4000 


DO 70 J = 1,NLOA 
PP eUCs ) ise MN 9) THEN 


MIN = UU(J) 


IND J 
COUNT = 1 
END IF 


CONTINUE 


P(e ecOuUNteaGeco ) THEN 


NNEG = I - 1 


GO TO 4000 
END IF 
V(I) = MIN 


VIND(I) = IND 


UU (IND) s2 0000 


CONTINUE 


CONTINUE 


C The Modified Depth First Search Algorithm. 


C Create all nodes out of the source and include them 


Gein a stack. 
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DO 80 I =70, 39 Gamat 


TIME(I) = 0 

SLACK (1) = 70 

LOAD(I) = 0 

PRED( 1) == 
80 CONTINUE 


DO. 90 1 = My 
XCOL(I) = 0.DO 
PATH(I-1) = 0 
CTBACK(I-1) = 0 
CTSON( f= 1) == 6 


20 CONTINUE 


It 
>) 


LENGTH 


( 
} 


LASTND 


TOP = 0 


CURED 


i 
© 


FROLD 


I 
o) 


TOLD = 0 


DO 100 I = NNEG,1,-1 


LOAD(NNEG - I + 2) = VIND(T) 
PRED(NNEG - I + 2) = 1 
STACK(NNEG - I + 1) = (NNEG - I + 2) 


TOPS) TOP stat 


LASTND = LASTND + 1 
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100 


e 


CONTINUE 


C Main loop to search for feasible schedules. 


@ 


3000 


tO 


12.0 


130 


140 


CURLD = STACK(TOP) 


IF ( 


LENGTH .EQ. NNEG ) THEN 
RCOST = 0.0D0 


DO 110 I = 1,LENGTH 


XCOL(LOAD(PATH(I))) = 1.0D0 
CONTINUE 
XCOEBANLOA + SHIP) @=51.0D0 


DO 120 I = 1,M 
Hee ce Oi saga ODO ) RCOSTe=—REOST + U(1) 
CONTINUE 
HR e@ot o-Gil 1 .0D—-4) RCOST = 0.0D0 
Pecos leno. ODO) )} THEN 
N = Nt+il 


BOs130°>F = 1,M 


A(I,N) = XCOL(T) 
CONTINUE 
XTIME(N) = TIME(PRED(CURLD) ) 


DO 140 J = 1,LENGTH 
SEQ(N,J) = LOAD(PATH(J) ) 
CONTINUE 


LSEQ(N) = LENGTH 
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RETURN 
END IF 
END IF 
TF CURED EQ. 0) SEnEN 
IF( MINRC -LT. 0.0D0) )eREIURT 
CTSHnir =-¢iohieaae 
SHIP = Shir ot 
EFASHIP .GT] NSH) SHEP ea 
IF( CTSHIP .EO. NSH)? feel 
WRITE (NOUT, 8000) 
TAU = TAU + 1 
WRITE(NOUT,8010) TAU 
210) At @ farsi 6 )16)16 
BND TF 
GO TO 2000 


END IF 


TOP = TOP. =) 


IF( N..GE. NN ) STOP "RUN CUT CheoEAecrE. 


IF( PRED(CURLD) .EQ. 1) THEN 
LASTND = LASTND - LENGTH 


LENGTH 


0 


END IF 
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i 0 


5000 


c 


IF( CTSON(LENGTH) .EQ. CTBACK(LENGTH) ) THEN 
DO 150 I = 1,LENGTH 
IF( PRED(CURLD) .EQ. PATH(I) ) THEN 


LASTND = LASTND - LENGTH + I 


LENGTH = I 
GO TO 5000 
END IF 
CONTINUE 
END IF 
PATH(LENGTH + 1) = CURLD 


C Compute the travel time to pick up another cargo. 


c 


c 


FROLD = LOAD(PRED(CURLD) ) 
TOLD = LOAD (CURLD) 
TT = 0 


CALL RTIME(NPOE,NPOD,MR,TR,IT,FROLD,TOLD,SHIP,NSH, 


&NLOA,TT,SPD) 


TIME(CURLD) = TIME(PRED(CURLD)) + TT 


C Verify if it is feasible to pick up another cargo. 


€ 
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170 


160 


ao, 


190 


IF ( TIME (CURDD) 32 LE CAC Spee 
CTBACK(LENGTH) = CTBACK(LENGTH) + 1 
LENGTH = LENGTH + 1 
CTSON (LENGTH) = 0 
DO 160 I = NNEG,1,-1 
DO 170 J = 1,LENGTH 
IF( VIND(I) .EQ. LOAD(PATH(J)) ) GO TO Wie 
CONTINUE 


LASTND = LASTND + 1 


LOAD (LASTND) = VIND(T) 
PRED(LASTND) = CURLD 

TOES -1OPr = ai 

CTSON (LENGTH) = CTSON(LENGTH) + 1 


CONTINUE 

DO 180 I = LASTND, (LASTND - CTSON(LENGTH) +i) 
STACK(TOP) = I 
TOE <> 70 OR Ga 

CONTINUE 

TOP = TOP + CTSON (LENGTH) 

ELSE 
LASTND = LASTND - i 
RCOST = 0.0DO 


DO 190 I = 1,LENGTH 


XCOL (LOAD( PATH( 2) ))) == sips 
CONTINUE 
XCOL(NLOA + SHIP) = 1.0D0 
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BOsz0e. 1 = 17M 
Petco seme O. ) RCOST = XCOST + U(I) 
200 CONTINUE 
MF(RCOST .GT. 9 -1.08-—4) RCOST = 0.0D0 
Pe kReCOST oon). soeomO .AND. LENGTH .GT. 
& INT(NLOA/NSH)-1 ) THEN 
IF ( CTBACK(LENGTH) .EQ. O ) THEN 
No = Not 1 


PO 210 19S eieM 


A(I,N) = XCOL(T) 
210 CONTINUE 
XTIME(N) = TIME(PRED(CURLD) ) 


DO 220 J = 1,LENGTH 


SEQ(N,J) = LOAD(PATH (J) ) 
BOM CONTINUE 
LSEQ(N) = LENGTH 


Pecos) wiih) MINRC  ) THEN 
MINRC = RCOST 
K = N 
ENDS IF 
END IF 


DO 230 I = 1,LENGTH 


XCOL(LOAD(PATH(I))) = 0.0DO 
XCOL(NLOA + SHIP) = 0.ODO 
230 CONTINUE 


ELSE 
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DO 240 I = 1,LENGTH 


XCOL (LOAD (PATH (I= Omer 
240 CONTINUE 
XCOL(NLOA + SHIP) = 0.0DO 
END IF 
CTBACK(LENGTH) = CTBACK(LENGTH) + 1 
GO TO 3000 
END iF 


IF (°N)2GE. LIMIT + aa eRe 
CTBACK(LENGTH) = 0 


GO 70.3000 


S000 FORMAT(//,6X,' TAU NOT FEASIBLE, SRNGREASE TAU) 
S010 "FORMAT 7, 6X7 NEW 2A0- = * | EAD 


END 


SUBROUTINE RTIME (NPOE,NPOD,MR,TR,IT, FROLD, TOLD, SHE 


&NSH,NLOA,TT,SPD) 


x This subroutine calculates travel times. 


IMPLICIT DOUBLE PRECISION( A-H,O-Z }, INTEGER (Sel siee 
PARAMETER( KK = 2, MM = 100, NN = 2000 ) 

DIMENSION SPD(MM) 

REA O71 5), TRS | 1S) 


INTEGER MR(100,KK) , TT, TOLD, SHIP EF Rerr 
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tl 
© 


pak 


C Calculating the travel time. 


IF(FROLD .EQ. 0) THEN 


TT = IDNINT((IT(SHIP,MR(TOLD,1)) + TR(MR(TOLD,1), 
& Mee OmD re) ey ten. SED (SHIP) }) 
ELSE 
TT = IDNINT((TR(MR(TOLD,1) ,MR(FROLD,2)) + 
& BoMiOoubD na, Mh | TOLD) 2) ))i7-(2ae* SPD(SHIP)) ) 
END IF 
RETURN 
END 


SUBROUTINE RDAYS (NLOA,NPOE,NPOD,NSH,TAU,TAUL, 


&TAUEST,SPD,TR,IT) 


* This subroutine calculates an initial estimate of the 


* number of days to complete the deployment. 


PARAMETER ( MM = 100 ) 
INTEGER NPOE,NPOD,TAU,TAUL, TAUEST,NLOA,NSH 
DAMENSTON TRiMoe15) , 20030,15) 


ReAG TRMIN, IT,7TR, lIMIN,SPMAX 


oy 


DOUBLE PRECMoLON ee pam 


COMMON /UNITS/ NIN,NOUT 


TAU = 1 

TAUL = 1 

TAUEST = 1 
C 


C Calculate the minimum distance to travel. 


C 
TRMIN “=S999o9999. 
DO 100 I=1,NPOE 
DO 200 29=15 NEOD 
IF(TR(I,J) .LT. TRMIN) TRMING= ere Go 
200 CONTINUE 


100 CONTINUE 


C Calculate the minimum travel distance from the initial 


C ships' locations to the POEs. 


C 
CEMEN = J9973990. 
DO 300 I = 1,NSH 
DO 400 J = 1,NPOE 
IF(IT(I,J) .LT. IPMIN) ITMIN — Ties 
400 CONTINUE 


300 CONTINUE 
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e 
C Calculate the maximum traveling ships' speeds. 
C 
SPEMAX = = 1. 
DO 600 J = 1,NSH 
IF(SPD(J) .GT. SPMAX) SPMAX = SPD(J) 


600 CONTINUE 


C Compute the average number of trips per ship. 


NTR = NLOA / REAL(NSH) 


C Calculate an estimate for TAU. 


Pee eee NTR see) 2) °* TRMIN) + ITMIN) / 
& (SPMAX*24.)) 
Covet ees. > + ( 0.8 * TAUL )) 


PAU S— TNieeeo. 7° * TAUL) + (0.3 * TAUEST) ) 


RETURN 


END 
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SUBROUTINE RESULT (JOUT,XB,U,C,A, TB (M,N Ober 


&SEQ,LSEQ, XTIME,NSH,NLOA) 


* This subroutine writes the solution to the output file 


IMPLICIT DOUBLE PRECISION(A-H, O-Z), INTEGER(I-N) 
PARAMETER ( MM =100, NN = 2000, ZERO = 0.0DO ) 
DIMENSION U(MM),C(NN) ,XB(MM) , IB (MM) 

INTEGER SEQ(NN,MM) ,LSEQ(NN) ,XTIME(NN) 

REAL A(MM,NN) 


COMMON 7UNDIs7 NIN, Nour 


PE CdOUT GE. 82) Co TOs s0 
WRITE (NOUT,8000) OBJ 
WRITE (NOUT 6005) “ETER 


WRITE (NOUT, 8010) 


G. Is-Xii} basic? 


c 
DO 30 I=1,N 
DO 10 J=1,M 
INDEX = J 
IF(IB(J) 2BO. I) Gomtomzs 
10 CONTINUE 


COnLo 30 


V2 


20 CONTINUE 
Welt etNeu) oOo) 1) xB 4 INDEX ) 
29, CONTINUE 
WRITE (NOUT, 8030) 
WRITE (NOUT,8040) (I,U(I) ,I=1,NLOA) 
WRITE (NOUT,8050) (I,U(I) ,I=(NLOA+1) ,M) 
DO 70 I=1,N 
DO 40 J=1,M 
TELE eho. ko CHEN 
Peer. Glee. D-2) AND. (1. .GT. 
& 2*NLOA+NSH) ) THEN 


DO 90 L=NLOA+1,M 


EPA A (ieee GT. -.9) 
& WRITE(NOUT,9010) (L-NLOA) 
90 CONTINUE 


WRITE (NOUT,9020) (SEQ(IB(J),K), 
& K = 1,LSEQ(IB(J)) ) 


WRITE(NOUT, 8090) IB(J) ,XTIME(IB(J)) 


END IF 
GOmTo. 60 
BN ar re 
40 CONTINUE 


C 
eux (l) is non basic. 


e 


i 


SUM = C(I) 
DO 50 J=1,M 
SUM = SUM + A(J,1) *U(J) 
50 CONTINUE 


GO 10770 


C 


C X({(I) is basic. 


€ 

60 CONTINUE 

70 CONTINUE 

80 CONTINUE 
IF (JOUT .£O. 2). WRITE(NOUT, c07¢) 
IF(JOUT .EQ. 3) WRITE(NOUT, 8080) 
RETURN 

e 


8000 FORMAT(//,6X,'OPTIMAL OBJECTIVE FUNCTION VALUE IS', 
SPs 5) 

8005 FORMAT(//,6X,'NUMBER OF ITERACTIONS = ',15) 

8010 FORMAT(//,17X,'OPTIMAL PRIMAL SOLUTION',/) 

8020. FORMAT (18x) 'X("). 13," )) =" busy 

8030 FORMAT(//,18X,'OPTIMAL DUAL SOLUTION',/) 

S040 FORMAT (18xX,"°U(" 712," )  =33F eae) 

8050 FORMAT(18X,'V(', 129°) =" fraaee) 

8070 FORMAT(//,6X,'PROBLEM IS UNBOUNDED FROM BELOW') 


8080 FORMAT(//,6X,'PROBLEM HAS NO FEASIBLE SOLUTION' ) 
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Pee CmmvOonrMmAl (6X, DURATION OF SCHEDULE °,14,' IS : ',I3,/) 

EneCme PORMAT(/,6X, SCHEDULE FOR SHIP ;: ',1I2,/,6X, ‘CARGOES 
&TO CARRY: ') 

9020 FORMAT (6X,14) 


END 


TO 


Or, 
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